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ABSTRACT 



This is the fourth paper in a series aimed at finding high-redshift quasars 
from five-color {u'g'r'i'z') imaging data taken along the Celestial Equator by the 
Sloan Digital Sky Survey (SDSS) during its commissioning phase. In this paper, 
we use the color-selected sample of 39 luminous high-redshift quasars presented 
in Paper III to derive the evolution of the quasar luminosity function over the 
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range of 3.6 < ;z < 5.0, and -27.5 < M1450 < -25.5 {Q = 1, Hq = 50 km s"^ 
Mpc"^). 

We use the selection function derived in Paper III to correct for sample 
incompleteness. The luminosity function is estimated using three different 
methods: (1) the estimator; (2) a maximum likelihood solution, assuming 
that the density of quasars depends exponentially on redshift and as a power 
law in luminosity, and (3) Lynden-Bell's non-parametric C~ estimator. All 
three methods give consistent results. The luminous quasar density decreases 
by a factor of ~ 6 from z = 3.5 to z = 5.0, consistent with the decline seen from 
several previous optical surveys at z < 4.5. The luminosity function follows 
ip{L) oc for z ~ 4 at the bright end, significantly flatter than the bright end 
luminosity function ipi^) oc L~^-^ found in previous studies for 2; < 3, suggesting 
that the shape of the quasar luminosity function evolves with redshift as well, 
and that the quasar evolution from z = 2 to 5 cannot be described as pure 
luminosity evolution. Possible selection biases and the effect of dust extinction 
on the redshift evolution of the quasar density are also discussed. 



Introduction 



Soon after quasars were first discovered, [Schmidt (1968)| found them to evolve strongly 
with redshift. In the following 25 years, a consensus has been reached that the number 
density of optical quasars peaks at 2; ~ 2.5 — 3, and declines towards both lower and higher 
redshifts. The luminosity function of optical quasars at z < 2.5 has been well studied 



through UV-excess and slitless spectroscopic surveys (e.g., [Boyle, Shanks fc Peterson 1988 



Hewett et al. 1991| , [Pei 1995| , |Boyle et al. 2000| ). At the high-redshift end, both multicolor 
surveys ( |Warren, Hewett fc Osmer 1994| , WHO hereafter; |Kennefick, Djorgovski fc de 
Uarvalho 1995| , DPOSS hereafter) and grism surveys ( [Schmidt, Schneider fc Gunn 199"B| , 



SSG here after) have shown that the number density of quasars declines rapidly from 2; ~ 3 
to 4.3. The sample sizes for these surveys are relatively small, so the rate of this decline and 
the shape of the high-redshift quasar luminosity function are not nearly as well-constrained 
as at low redshift. The quasar luminosity function provides important constraints for 
models of quasar evolution (e.g., [Efstathiou fc Rees 1988] , [Turner 1991| , [Haiman fc Loeb 
1998| , [Haiman fc Hui 2000| ), and the nature of the high-redshift UV ionizing field ( [Madau 



Haardt fc Rees 1999[) . 

In Paper III ( [Fan et al. 2000a[) , we presented a sample of 39 quasars at 2 > 3.6 and 
i* ^ 20, selected from 182 deg^ of Fall Equatorial Stripe multicolor imaging data taken 
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during the commissioning phase of the Sloan Digital Sky Survey (SDSS, [York et al. 20001 , 
see also [bukugita et al. 199^ , |Gunn et al. 199^ ). The selection completeness function of this 
multicolor sample is calculated based on model quasar colors. In this paper, we derive the 
high-redshift quasar luminosity function based on the quasar sample and selection function 
presented in Paper III. In §2, we use three different statistical approaches to find the 
luminosity function: a 1/^4 estimator; a maximum likelihood solution; and Lynden-Bell's 
non-parametric C~ estimator ( [Lynden-Bell 1971] ). In §3, we discuss the rate of decline of 
the quasar spatial density towards high redshift and the shape of the quasar luminosity 
function based on the results of §2, and compare these new results with those of previous 
studies. In §4, we discuss possible selection biases on the evolution of the luminosity 
function and theoretical implications of our results. 



2. Derivation of Luminosity Function 

The differential quasar luminosity function, z), is defined as the number of quasars 
per unit comoving volume, per unit luminosity as a function of luminosity and redshift. It 
can also be written as \1/(M, z), the number of quasars per unit comoving volume, per unit 
absolute magnitude: 

■ dL 



^(M, z) = ^{L,z) 



dM 



0.92L^(L, z). 



Throughout the paper, we use 2;) for the mathematical derivations, while we present 
our final results in terms of the more familiar \1/(M, z). The quasar luminosity function can 
also be presented in the cumulative form: z) = "^{L', z)dL'. 

In general, the redshift and luminosity dependences of the quasar luminosity function 
are not separable. In limiting cases, z) can be represented either in the form of pure 
density evolution: 

^mz)=tlj{L)p{z), (2) 

where the luminosity function remains the same shape, but with different normalization at 
different redshift; or in the form of pure luminosity evolution 

^mz) = nL/giz)]/giz), (3) 

where g{z), the characteristic luminosity, is a function of redshift while the normalization 
of the luminosity function is a constant. It might be better presented by a combination of 
the two cases above (e.g., |Koo &: Kron 1988]) : 



<f{L,z)=piz)nL/g{z)]/g{z). 



(4) 
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For -2^3, pure luminosity evolution provides a good fit to the quasar luminosity 
function (e.g. ^chmidt fc Green 1978| , [Boyle et al. 2000| ). At 2; > 3, existing quasar 



samples are still too small to distinguish these models. In general, the evolution can be 
more complicated. For example, the luminosity function at low and high redshifts may be 
determined by very different physics and have totally different shapes. Note that if the 
luminosity function is a power law, \1'(L, 2;) oc L^, where f3 does not depend on z, there 
is no characteristic luminosity, and pure luminosity and pure density evolution cannot be 
distinguished. 

The determination of the high-redshift quasar luminosity function from a flux-limited 
sample faces several difficulties: (1) the selection function of a multicolor survey is a strong 
function of luminosity, redshift and the shape of the quasar spectral energy distribution 
(SED, see Paper 111); the selection bias has to be corrected for carefully. (2) The 
high-redshift quasar samples are typically small (the largest complete sample of quasars 
at 2; > 4 prior to this work has 10 objects). The small sample size makes it difficult to 
parameterize the luminosity function, and the result is also strongly affected by the choice 
of binning. (3) As the redshift increases, an optical survey samples different parts of the 
quasar's intrinsic spectrum. When comparing results from low and high redshift samples, 
the K-correction is large and quite uncertain, due to the lack of statistics on the continuum 
shape between the rest-frame UV and optical region of high-redshift quasars. 

The luminosity function can be derived using either non-parametric or parametric 
methods. The most commonly used non-parametric method is the 1/Va estimator (e.g. 
Avni fc Bahcall 1980| ). We present the 1/Va results in §2.1. In §2.2, we generalize the test 



of correlation between luminosity and redshift distributions described in pUfron fc Petrosian 



1992)1 and [Maloney fc Petrosian (1999)| , by including the selection function. We show that 



for the limited range of redshift and luminosity our survey covers, the two distributions are 
not correlated. Thus in §2.3, we assume that the luminosity and redshift dependences are 
separable, and that the quasar density follows a power law in luminosity and depends on 
redshift exponentially. The maximum likelihood solutions are derived following the method 
of [Marshall (1985)| , where each quasar is treated as a 5 function in parameter space, and 



the results are not subject to binning. When the density and luminosity evolution can be 
separated, one can also use the Lynden-Bell's C~ estimator ( [Lynden-Bell 1971 ) to calculate 



the non-parametric marginal distributions of the cumulative luminosity function along the 
luminosity and redshift directions. In Appendix A, we show how to generalize the C~ 
estimator to include the survey selection function. The luminosity function estimated using 
the C~ estimator is presented in §2.4. 



The results throughout this paper are presented for two cosmological models: (1) a 
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model with = 1 and Hq = 50 km Mpc~^, which we refer to as the f2 = 1 model. Most 
of the previous studies of the quasar luminosity function have used this cosmology. (2) a 
A-dominated flat model with Q = 0.35, A=0.65, and Hq = 65 km s~^ Mpc~^ ( pstrikef] 
fe Steinhardt 1993 , Krauss fc 'I'urner 1995 , [Bahcall et al. 11?^ ), which we refer to as the 



A-model. The luminosity function is presented in term of M1450, the absolute AB magnitude 
of the quasar continuum in the rest-frame at 1450 A. Assuming a power law /j, oc u", 
Schmidt et al. (1995) find: 

Mb = M1450 + 2.5a log(4400/1450) + 0.12, (5) 

where the effective wavelength of the Kron-Cousins B band is 4400 A, and the factor 0.12 
comes from the zero point difference between the AB and Vega-based magnitude systems 
for quasar-like spectra. 



2.1. 1/Va estimate 

We first calculate the space density of quasars in different luminosity and redshift bins 
using the 1/Va method following the discussions in WHO and DPOSS. The available volume 
Va is defined for a quasar of given luminosity and SED as the comoving volume over which 
the quasar could have been detected by the survey. In the presence of a selection function, 
the available volume is weighted by the selection probability. If one is to determine the 
luminosity in luminosity and redshift bins Az, AL, the available volume for a quasar with 
given L, z and SED is: 

K = / p{L,z,SED)^dz, (6) 
JAz dz 

where p{L, z, SED) is the selection probability of the quasar as a function of redshift, 
luminosity and the SED. The luminosity function and its statistical uncertainty can be 
estimated as: 



a 



? WAl) 



1/2 



(8) 



where (L) and (z) are the average luminosity and redshift over the bin, the sum is over all 
the objects in the bin, and the available volume is calculated for each object. 

In Paper III, two flavors of selection probability are calculated: (1) the probability for 
a given quasar SED shape, p[L, z, a, EW (Lya+NV)], where a is the slope of the power law 
continuum and EW(Lya;+NV) is the rest-frame equivalent width of the Lya+NV emission 
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line. These two parameters characterize the shape of the SED. The selection probabihty 
of each quasar in the sample is listed in Table 4 of Paper III. (2) The average selection 
probability p{L,z) (§5.3 in Paper III), averaged over the intrinsic distribution of a and 
EW. In Paper III, we assume that this distribution is the product of Gaussians in a and 
EW(Lya + NV), and find that their mean and standard deviations are: a = —0.79 ± 0.34 
and EW(LyQ; + NV) = 69.0 ± 18.3 A. In Figure 1, we present the 1/Va results of the 
quasar luminosity function using both flavors of selection probability under the — 1 
and the A-model cosmologies. Figure 1(a) and (b) are for the Q — 1 model. In Figure 
1(a), the density is calculated using the average selection function p{L, z), while in Figure 
1(b), the density is calculated using the selection probabihty of each quasar for its SED 
type, p[L, 2;, a, EW(LyQ;+NV)]. Similarly, Figure 1(c) and (d) are for the A-model. In each 
case, we divide the sample of 39 quasars in 3 redshift bins: 3.6 < z < 3.9 (18 quasars), 
3.9 < z < 4.4 (14 quasars) and 4.4 < z < 5.0 (7 quasars). The redshift bins are further 
divided into two or three luminosity bins. There are 3-6 quasars in each bin. 

The different methods of treating the selection function give consistent results. 
This can be understood as follows. For simplicity, we ignore the dependence on EW, 
z and L, and assume that the quasar distribution only depends on a. For an intrinsic 
distribution function f{a) and a selection function p{a), the observed distribution will be 
n{a) = p{a)f{a). In the discrete case, n{a) = J2i ^{<^ — o^i)- The effective density of quasars 
can be estimated as: 

Peff = ^ [ f{a) da 

1 f p{a)f{a) da 



Va J p{a) 
1 f n(a) da 



K i Pia 

Alternatively, if we define the average selection probability: 



Va J p{a) 



/p(a)/(«) da 
^ ina)da ' ^'^^ 

then the density can also be estimated as: 

fp(a)f(a)da N 
Pes = =77 = =77, (11) 

pVa pVa 

where N is the total number of objects observed. Therefore, if, as we have assumed, 
the SED shape is not a function of redshift or luminosity, the two approaches only differ 
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in the different weightings apphed to the observed n{a) distribution; in Eq. (11), each 
object is given the same weight and normahzed by the average selection probabihty. The 
two methods should give consistent results if the sample is large enough and all the SED 
shapes are properly sampled. However, as we showed in Paper III, there are regions in the 
(L, z,a,EW) parameter space poorly sampled by the selection criteria {p <C 1). In this 
case, the average probability p{L, z) is much better behaved, and we will use it for most of 
the analyses which follow. 

Due to the small number of quasars in each bin, the appearance of Figure 1 changes 
somewhat with different binning, but the general trends are obvious: (1) the overall density 
of quasars declines rapidly with redshift. For the = 1 model, the density of quasars at 
-^1450 ~ —26.5 decreases from ~ 2 x 10"*^ Mpc~^ mag"^ aX z = 3.75 to ~ 6 x 10~^ Mpc~^ 
mag~^ at z = 4.75. (2) the luminosity function rises towards lower luminosity by a factor of 
~ 4 per magnitude. These results are very similar for the A-model. At z ~ 4, the behavior 
of the n = 1 and A models are quite similar, as VL{z) 1 when the redshift becomes large. 

2.2. Test of Correlation between the Luminosity and Redshift Distributions 

In this subsection, we test whether the intrinsic luminosity and redshift distributions of 
the quasar sample are correlated or whether they can be treated separately (as in the case 
of pure density evolution, Eq. 2). At low redshift, the luminosity function can be fit with a 
double power law (e.g., [Boyle, Shanks fc Peterson 1988|) . The sample in this paper covers a 
limited range of absolute magnitude (—27.5 < M1450 < —25.5) and a limited redshift range 
(3.6 < z < 5.0). If, in this case, the luminosity function can be fit with a single power 
law whose slope is independent of the redshift, the luminosity and redshift distributions of 
our sample are not expected to be correlated. This correlation test is also crucial for the 
non-parametric C~ estimator in §2.4, which requires that the data be expressed in terms of 
two uncorrelated variables. 

Correlation tests for truncated data sets (such as fiux-limited samples) have been 
developed by Efron & Petrosian (1992) and [Maloney fc Petrosian (1999)| . In our survey, 
there is essentially no bright-end truncation in the selection (the saturation limit is i* ~ 14, 
four magnitudes brighter than the brightest quasar in the sample). Thus, our data set 
{Lj, Zi} is one-side truncated in luminosity. We first consider the case in which the selection 
function is a step function and the selection boundary is sharp: 




(12) 



- 9 - 



where L~{z) is the hmiting luminosity at a given redshift. In this case, following the 
discussion in |Maloney fc Petrosian (1999)| , one first defines the comparable or associated set 
for each object i: 

J, = {j:Lj>Li,L-<L,} (13) 

consisting of all objects which are brighter than the object in question and which would 
be selected by the survey if they had the same luminosity as the object in question. This 
is the same set that will be used in the C~ estimator below. It is the luminosity and 
volume limited subset of the sample that can be constructed for each object. If z and L are 
independent, the rank i?, of Zi in the comparable set, defined as 

Ri = number of j, such that {j G Jf, Zj < Zi}, (14) 

should be distributed uniformly between and Ni, where Ni is the number of points in 
the comparable set for each object (not including the object in question). The expectation 
value of Ri, is Ei = {l/2)Ni and the variance Vi = {1/12)N!^. Therefore, one can construct 
Kendell's r statistic: 

r-^. ,15) 

Foi' I ''"I ^ 1) the luminosity and redshift would not be correlated at ~ la level and could be 
treated independently. 

However, the discussion above only applies to a sample with a sharp boundary. For 
the quasar survey we discuss in this paper, even though the selection criteria have a sharp 
cut at the limiting apparent magnitude {i* < 20.05 for gri selected quasars and i* < 20.2 
for riz selected quasars, see Paper III), the selection function smoothly approaches zero 
at the limiting luminosity (Figure 8 of Paper III). There is no one-to-one correspondence 
between the apparent i* magnitude and the luminosity of the quasar at a given redshift, 
due to scatter in the quasar SED shape and photometric errors in the SDSS measurements. 
We thus generalize the correlation test to include the selection function p{L,z). First, we 
define a generalized comparable set: 

J, = {j : > L,}, (16) 

which includes all the objects more luminous than the object in question, since all might 
have a selection probability between and 1. We then define the total number in the set 
by weighting each point j in the comparable set Jj. This weight is proportional to p{Li, zj), 
the selection probability if the object j had the same luminosity Li as the object in question 
(object i), and it is proportional to the inverse of its own selection probability p{Lj, zj). It 
is useful to define a quantity T^: 
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The definition of the rank changes accordingly: 

^_^p^^ if^.<^.. (18) 

The expectation value of the distribution Ri is Ei = (l/2)Tj, and the variance 

Vi = (l/12)(Tj)2. The test statistic r is defined in the same way as in Eq. (15). Note that 

these definitions return to its normal form if the selection function has a sharp boundary, 

since p{Li,Zj) — when Li < LJ . In the case where the selection boundary is smooth, 
the ranking is weighted by p[Li, Zj)/p{Lj. Zj). The selection probability always increases 
with larger luminosity and smaller photometric error. Therefore, since p{Li, Zj) < p{Lj,Zj) 
for Li > Lj, the weight is always between and 1. In Appendix A, we prove that the 
comparable set defined in this way gives the correct estimator for the luminosity 
function in the presence of a selection function. 

We calculate the ranking statistic r of the SDSS Equatorial Stripe quasar sample using 
the average selection function in Paper 111. For the Q = 1 model, r = —0.20; for the A 
model, r = —0.15. Therefore, for the limited redshift and luminosity range this sample 
covers, the two distributions can be treated independently. For the full quasar sample 
of the main SDSS survey which will cover a much larger redshift and luminosity range, 
the possibility of luminosity evolution will need to be introduced. The ranking technique 
presented here can then be used to determine the functional form of the luminosity evolution 
{g{z) in Eq. 3) of the quasar population (e.g. Maloney & Petrosian 1999). 



2.3. Mciximum Likelihood Fits 

In this subsection, we model the luminosity function using the result above that the 
redshift and luminosity dependences are separable. Imagine dividing up the sample into 
small bins in L and z, Li and zj, with sizes AL and Az. The likelihood that riij quasars are 
found in the (Lj, zj) bin can be written as : 

where /lij is the average number of quasars expected in the bin: 

r r dV 
IXij= ^{L,z)p{L,z)—dL dz, (20) 

J Jal,Az az 

where ^{L,z) is the differential luminosity function, and the integral is over the bin. The 
maximum likelihood solution is obtained by minimizing the function 

5 = -21n£, (21) 
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with respect to the parameters describing the luminosity function. 

One could find the maximum likehood solution by fitting the 1/Ki results in §2.1, but 
the result would be strongly affected by the choice of binning for a small sample. However, 
in the limit of an infinitesimal bin size, riij is either and 1, and the result is not dependent 
on binning. Following [Marshall (1985)| , the likelihood function can then be written as: 

^ r r dV 

S = -2J2 H'^iLi, z,)p{Li, z^)] + 2 J J ^(L, z)p{L, z)^dL dz, (22) 

i 

where the sum is over all quasars in the sample. The integral is over the entire (L, z) range 
of the survey, and is equals to the number of quasars expected in the survey for a given 
luminosity function. This term provides the normalization for the likelihood function. 

The likelihood function can be expressed in a slightly different way (see also §5.2 in 
Paper III): 

\^JJ^{L,z)p{M,z)%dLdz ^ ' 



S' = -2 In/:' = -2Y,H'^ {Li, Zi)p{Mi,Zi)] + 2N In (^JJ "^{L, z)p{L, z)^dL dz^ 



(24) 



where N is the sample size. As in Eq. (22), the integral term provides a normalization. 
These two expressions are closely related. If we define N' = JJ \I'(L, z)p(L, z)^dL dz, then 
for any parameter x, it can be shown that 

?^^?^^2?f(^-l). (25) 



dx dx dx \N' 

Thus if we constrain the total number of objects expected based on the maximum likelihood 
solution to be equal to the total number actually observed, A^' = A^, Eqs. (22) and (24) are 
maximized at the same parameters. In practice, they give almost identical results. 

The calculation in §2.2 shows that the luminosity and redshift distributions of the 
quasar sample are not correlated. We assume that the quasar density is a power law function 
of the quasar luminosity : ip{L) oc L^, or, expressed in terms of absolute magnitude, 
■ip{M) oc iQ-0-4:(P+i)M _ Following SSG, we further assume that the density at a given 
luminosity declines exponentially with redshift p{z) oc 10""^^. SSG express their results in 
the form of a cumulative luminosity function 

log$(z, < Mi45o) = A-B{z-3) + C(Mi45o + 26), (26) 

This is equivalent to the differential luminosity function 

•qj* 

"^{Z, Mi450) = io0.4[Mi450+26-a(.-3)](/3+l)' (^7) 
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where ^* = io-o-4inio(/3+iM^ p ^ _2 _ 1, and a = -2.5B/{P + 1). 

As discussed in §5.2 of Paper III, the average probabihty p(M, z) depends on the 
underlying distribution of quasar SED parameters a and EW(Lya+NV). The hkehhood 
function (Eqs. 22 and 24) has to be maximized not only with regard to the luminosity 
function, but also to the intrinsic distribution of a and EW(Lya+NV) at the same time. 
In practice, however, the maximum likelihood solutions of the {a, EW) distribution and 
of the luminosity function are only weakly correlated. The solutions can be found by first 
assuming the a and EW distributions based on their weighted average {albeit biased by 
selection effects) to derive the first-order luminosity function and then iterating the results. 

Using the likelihood function in Eq. (22), we solve for the parameters of the cumulative 
luminosity function (Eq. 26). For the = 1 model, 

log$(2,< M1450) = (-7.24 ± 0.19) - (0.48 ±0.15)(2- 3) + (0.63 ±0.10)(Mi45o + 26), (28) 

where $ is in units of Mpc~^. From the maximum likelihood fit, we also find a = — 0.79±0.34 
and EW(Lya+NV) = 69.3 ± 18.0 A, as reported in Paper III. The best-fit parameters 
in both cosmological models and in both cumulative and differential forms are given in 
Table 1, and these models are plotted as dashed lines on Figures 1 and 2. The error bars 
on the parameters are estimated both by finding the parameters that yield S = Smin + 1 
(which gives the correct la errors for a Gaussian distribution), and by using the bootstrap 
technique. The results from the two methods are consistent with each other within 10%; 
the values given in Table 1 are from the Smin + 1 method. We find that the best-fit 
luminosity function predicts a total number of // "^{M, z)p{M, z)^dMdz = 38.7(38.6) 
quasars observed in the survey area for the Q = 1 (A) model, compared to the 39 quasars 
actually observed in the survey. 

2.4. Non-Parametric Determination : Lynden-Bell's C Estimator 

If the luminosity and redshift distributions of objects in a sample are uncorrelated 
(§2.2), the bivariate luminosity function \E'(L, z) can be separated: \E'(L,2;) = ilj{L)p{z), 
where ip{L) and p{z) are the marginal distributions of the luminosity function in the 
redshift and luminosity directions. In this case, the cumulative marginal distribution in the 
luminosity direction, (f){L) = J ijj{L)dL can be estimated using Lynden-Bell's C~ estimator 
( ILynden-Bell 1971|) : 

^{Lj)=<P{L^)f[{l + l/N,), (29) 

k=2 
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where A^^ is the number of objects in the comparable set of the object k as defined in Eq. (13), 
and the objects are sorted according to their luminosities: Li > ... > > Li > ... > L^- 
Similarly, the cumulative marginal distribution in redshift a{z) = J p{z)dz can be written 
as: 

i 

a{zi) = a{z^)X{{l + l/Mk), (30) 

k=2 

where is the number of objects in the box L > L^. and z < z^axiLk), and ZmaxiL^) is 
the maximum redshift the survey can detect at the luminosity Lk- Therefore, Mk is the 
size of the comparable set in the redshift direction, defined in a similar fashion to Nk. Here 
the objects are sorted in the order of their redshifts: Zi < ... < Zi_i < Zi < ... < z^. If 
the L and z distributions are correlated, one should find a set of parameters x = x{L, z) 
and y = y{L,z), whose distributions are uncorrelated, and use the C~ estimator to obtain 
the marginal distributions in x and y (e.g., Maloney & Petrosian 1999). Note that these 
equations only give the shape of the distribution, but not the overall normalization of the 
function. The distributions can be normalized by requiring tha the total predicted number 
of objects equal that observed: 

/•OO rOO 

Nobs= ilj{L)a[zmax{L)]dL = p{z)(t)[Lmin{z)]dz (31) 
Jo Jo 



As shown in [Petrosian (1992)1 , a variety of non-parametric methods can be reduced to 



Lynden-Bell's C method ( [Lynden-Bell 1971]) in the limiting case of one object per bin. 



The C~ estimator makes efficient use of the data by deriving the marginal distribution for 
the uncorrelated variables, and does not require binning when calculating the cumulative 
distribution. Obviously, binning is still needed when deriving the differential distribution. 
But in this case, rather than binning the data in both L and z as for the 1/Va estimator, 
we only need to bin the data one axis at a time. 

For survey with a complicated selection function p{L,z), the definition of the 
comparable set should be changed (Appendix A). For the luminosity direction: 

A^j = ^ ^~7j~^—^, where the sum extends over Lj > Li, (32) 
and for the redshift axis: 



Mi = y^ ^-tt~~Tj where the sum extends over Zk < Zi. (33) 
Y P{Lk, Zk) 

With these definitions, Eqs. (29) and (30) remain valid. 

We use the equations above to derive the marginal distributions for the SDSS sample. 
For the shape of the luminosity function, we calculate the quantity ilj{M)/(j){M < —25.5). 
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the number density of quasars as function of redshift, normalized by the total number 
at M < —25.5. For the redshift evolution, wc calculate p{z, M < —25.5), the total 
spatial density of quasars with Afi45o < —25.5 as a function of redshift normalized via 
Eq. (31). Note that p{z,M < -25.5) = p{z)(f){M < -25.5), therefore the product of the 
two quantities gives the differential luminosity function: these distributions are shown in 
Figure 2 for both the ^2 = 1 model and the A model, and compared with the maximum 
hkelihood results. The data are binned to have roughly the same number of objects per 
bin. Prom Figure 2, we find ^/j{Mu5o) oc iq-o-^sa^wso^ and p{z) oc IQ-^-^^^ for the Q = 1 
model. These two coefficients are —0.63 ± 0.10 and —0.48 ± 0.15 in the maximum likelihood 
fits (the dashed hues in Figure 2). Comparing Figures 1 and 2, the C" results benefit 
from the smaller bin size and have higher signal-to-noise ratio because they use marginal 
distributions. 



3. Evolution of the Quasar Luminosity Function at High Redshift 

In §2, we calculated the quasar luminosity function from the SDSS Equatorial Stripe 
sample using both parametric and non-parametric methods. This quasar sample includes 

39 quasars at i* <^ 20 and covers 182 deg^. Among them, 18 objects are at 2; > 4.0 and 
six at z > 4.5. This is the largest complete quasar sample at z > 3.6 to date. We derive 
the luminosity function over the range 3.6 < z < 5.0 and —27.5 < M1450 < —25.5 (for the 
Q = 1 model). The three methods (the 1/Va estimator, the maximum likelihood solution 
and Lynden-Bell's estimator) yield consistent results: 

1. Over the redshift and luminosity range considered, the distribution of z and L are 
uncorrelated, therefore, the luminosity function can be separated: ^{L, z) = ijj{L)p{z). 

2. The quasar number density declines rapidly towards high redshift. It can be fitted by 
an exponential decline: p{z) ~ iq-0-5z ^ g-i.i5z. ^j^^ spatial density drops by a factor 
of ~ 3 per unit redshift. 

3. The quasar number density rises towards fainter luminosity as ijj(M) ~ 10°'^M, or 
ip{L) ~ L~^-^. The density increases by a factor of ~ 4 per magnitude. 

In this section, we first compare these results with previous samples at high redshift, 
and then discuss the shape of the high-redshift luminosity function. 
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3.1. Comparison with Other Surveys at z > 3.6 

The quasar luminosity function a.t z > 3.6 has been calculated for several previous 
optical quasar surveys (WHO, SSG, DPOSS). The selection of all these surveys is based 
upon the observations in the rest-frame UV, as in our survey. However, some of the results 
are expressed in the absolute i?-band magnitude Mb in the rest-frame. The UV magnitudes 
are converted to B-band magnitudes assuming a power law continuum fi^ oc z/~°-^ (SSG), 
which we adopt for consistency. From Eq. (5), we have Mb = M1450 — 0.48. We re-write 
our maximum likelihood results in the f2 = 1 model: 

log<l>(z,< Mb) = (-6.91 ±0.19) - (0.48 ± 0.15)(2 - 3) + (0.63 ± 0.10)(Mb + 26) (34) 

where $ has units of Mpc^'^. Since all the surveys to which we will compare are selected 
from the rest-frame UV fluxes at similar redshifts, uncertainties in the K-correction will not 
affect the comparison (see also §4.1). 

SSG derive the quasar luminosity function in the redshift range 2.75 < z < 4.75 for 90 
quasars at r ^ 21, selected by their Lya emission in the Palomar Transit Grism Survey. 
The survey covers 61 deg^, and includes 20 quasars at z > 3.6 and 9 quasars at z > 4.0. For 
the = 1 model, SSG find: 

log$(z, < Mb) = -6.84 - 0A3{z - 3) + 0.75(Mb + 26). (35) 

In Figure 3, we compare the quasar evolution from the SDSS and from the WHO and 
SSG surveys. We also show the low-redshift results from the 2dF survey ( |Boyle et al. 2000| ). 
The new SDSS result is consistent with previous surveys. It shows very similar redshift 
evolution to that of SSG, and this trend continues towards z ~ 5. The SDSS luminosity 
function is somewhat flatter than that of SSG, although the difference is only at the 1-cr 
level, and is based on rather small sample sizes. Note, however, that the median redshift of 
the SDSS sample {z ~ 4) is considerably higher than that of SSG's (2 ~ 3.3), which could 
have an effect on the luminosity function slope (see §3.2). 

The consistency of these two surveys is important because they used totally different 
selection methods and have very different selection functions: SSG is a slitless spectroscopic 
survey whose selection is based on the detection of emission lines in grism spectra. It is 
therefore is biased towards quasars with strong emission lines, but is not selected on the 
continuum shape of the quasars. On the other hand, the SDSS sample selection is based on 
broad-band colors. It depends on the strength of the emission lines rather weakly, but is 
biased towards objects with stronger continuum breaks and/or bluer continua, and is not 
sensitive to red quasars (a ^ —1.6, see the discussion in Paper III). 
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The WHO survey (|Warren, Hewett fc Osmer 1994 ) and the DPOSS survey ( KenneticF 



Djorgovski fc de Carvalho 1995| ) used broad-band colors to select quasar candidates. These 



studies are based on photographic photometry from Schmidt plates, which typically have 
photometric errors ^ 0.10 mag. They use a very similar simulation technique to Paper III 
to correct for selection effects, although they suffer from a larger selection incompleteness 
(~ 50%). The WHO survey covers an effective area of 43 deg^ down to rrior = 20. It 
contains 86 quasars at z > 2.2, including 8 at 2; > 3.5 and 2 at 2; > 4. They found that the 
quasar number density drops by a factor of ~ 6 from z = 3.3 to 2; = 4.0, although given 
the small number of objects in their sample, the rate of evolution is poorly constrained. 
The DPOSS survey includes 10 quasars at z > 4 and 16.5 < r < 19.6, covering 681 deg^. 
They derive the quasar luminosity function at their median redshift z = 4.35. In Figure 4, 
we compare the cumulative quasar luminosity function from the SDSS sample at 2; ~ 4.3 
with the results from these three previous surveys. The results are all consistent within 
the error bars. Note that although the two previous multicolor surveys needed a much 
larger selection correction than the SDSS sample, the final corrected results agree with one 
another quite well. 

The comparisons above show that the decline of the number density of luminous 
quasars from optical surveys is unlikely to be strongly affected by either the selection 
technique or the correction for the selection functions. In §4.1, we discuss the effect of 
K-correction and extinction on the evolution of high-redshift quasar densities. 



3.2. The Shape of the High-redshift Quasar Luminosity Function 

Using the SDSS sample, we find that at 2 ~ 4 and Mb < —26 , the shape of the quasar 
luminosity function can be fitted with a power law 

^{L) oc L-2-5±o-25, (36) 

where the power law index is the average of the results from maximum likelihood and C~ 
estimators. 

At 2; < 3, the quasar luminosity function is often fitted with a double power law: 

"^{Mb, z) = ^Q0.4[(/3i+l)(Affl-Af*(^))] ^ iQOAm+l)iMB^M*g{z))] ' (^7) 

In the expression above, pure luminosity evolution is assumed, with a redshift dependence 
given by the evolution of the characteristic luminosity M'^{z). Boyle et al. (2000) find that 
for the 2dF sample, the normalization of the luminosity function '^*{Mb) = 1.1 x 10~^ 
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Mpc ^ mag the faint end slope Pi — —1.58, and the bright end slope P2 — —3.43 (the 
Q—1 model). The error bar on the slopes Pi and P2 are ^0.1. 

If we assume that the luminosity function in Eq. (37) applies to 2; ~ 4, i.e., the shape 
of the luminosity function remains the same while only the characteristic luminosity M^{z) 
varies, we find that in order to match our results at Mb ~ —26, we must have ~ —24.7, 
more than one magnitude fainter than our survey limit. Therefore, the SDSS results in this 
paper clearly probe the bright end of the quasar luminosity function. 

The results on the bright-end slope of the low-redshift quasar luminosity function from 
various surveys appear to be quite robust: (f){L) oc L^, with P ~ —3.4 — 3.5 (Boyle et al. 
1988, Hewett et al. 1991, Pei 1995, Boyle et al. 2000), much steeper than our high-redshift 
results {P ~ —2.5). At low redshift, the quasar number density drops by a factor of ~ 10 
per magnitude at the bright end, while at ^ ~ 4, it only drops by a factor of 4. 

We can further test the significance of this discrepancy. We derive a new maximum 
likelihood solution from the SDSS sample by forcing the slope of the luminosity P = —3.43 
but allowing the normalization \['o and the redshift evolution a to vary in the luminosity 
function (Eq. 27). This is shown as the dashed line in Figure 5. The luminosity function 
from the SDSS sample using both the C~ estimator and the maximum likelihood estimator 
from Figure 2 are shown as the points and solid line in the figure. It is evident that this 
steep slope is not consistent with the data. The statistic S (Eq. 22) of the maximum 
hkelihood solution follows a distribution. In this case, the hkelihood ratio for a specific 
Po can be estimated by calculating AS — S{Po) — Smin, where Smin is the minimum S value 
when all three parameters are allowed to vary, and S{Po) is the S value if P is fixed to be 
Po and the other two parameters are allowed to vary. We find that by forcing P = —3.43, 
AS = 8.9 with two degrees of freedom. This corresponds to a probability p = 1%, so the 
two slopes are inconsistent at ~ 2.5-0" level. 

We can use another method to estimate the significance of the inconsistency between 
low and high-redshift slopes: there are six quasars at M1450 < —27 in the sample. Using 
the model with P = —3.43, only ~ 1.7 quasars are predicted over the whole survey. The 
probability of > 6 quasars being observed in our survey is only 0.8% from Possion statistics. 

It is unlikely that this fiattening arises because we are missing a large fraction of low 
luminosity quasars. The photometric errors are quite small and the selection function is 
very uniform as a function of magnitude (Figure 8 of Paper III) for all objects a few tenths 
of a magnitude brighter than the survey limit. Figure 5 shows that the slope is much fiatter 
than P — —3.43 even over the three most luminous magnitude bins. Thus the bright end 
slope of the quasar luminosity function is considerably flatter at 2; ~ 4 than at 2; ~ 2. 
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Note that SSG find il^{L) oc L~^-^, while their median redshift is somewhat lower. The 
fiattening of the bright end luminosity function at high- redshift, if confirmed using a larger 
sample, indicates that the evolution of quasar luminosity function at z > 3 cannot be pure 
luminosity evolution, unlike the case for z < 3. Instead, the shape of the quasar luminosity 
function evolves with redshift as well. 



4. Discussion 

Figure 3 shows that the comoving quasar number density peaks at z ~ 2.5, and drops 
by a factor of ~ 20 from z ~ 2.5 to 5. This decline is evident in all optical high-redshift 
quasar surveys to date. However, these results might be biased by three factors: the survey 
incompleteness; uncertainties in the K-correction when comparing low and high redshift 
results; and the possibility that a large fraction of quasars are not detectable in optical 
survey due to dust extinction. 

In Paper III, we show that the survey selection criteria are only sensitive to a range 
of the quasar SED shapes, and that the incompleteness is a function of redshift. When 
calculating the luminosity function, we first correct this incompleteness by assuming that 
the distribution of the continuum power-law index o; is a Gaussian with mean —0.79 and 
standard deviation 0.34. At 4.5 < z < 5.2, the survey is sensitive for quasars even with 
a < —2.0, although no quasars with a < —1.6 are detected in the sample. At z < 4.5, the 
survey is not sensitive to quasars with a ^ —1.6. Therefore, the density at ^ ~ 4 could be 
underestimated if there was a significant contribution from quasars at a ~ —2.0 that we 
failed to detect. However, since we are not biased against these red quasars at the highest 
redshift bin {z > 4.5), the strong decline in the quasar density seen at 2; ~ 5 will not be 
affected by them. Since our selection is not sensitive to very red quasars {a < —2.5), if 
there were a large population of such objects only at high redshift, the inferred density 
decline at high redshift from our color-selected sample might be an overestimate. But as 
we pointed out in §3.1, the same decline is observed in a spectroscopically selected sample 
(SSG), which argues that neither survey suffers from large incompleteness. 

Low-redshift quasar luminosity functions usually are expressed in the rest-frame B 
band magnitude Mb- When comparing with high-redshift results, which are measured in 
the rest-frame UV (e.g. M1450), a K-correction is needed. It is customary to assume a 
a — —0.5 power law continuum for the K-correction. An incorrect K-correction might result 

in a large difference in density. In fact, we find that the UV continua of quasars at ^ ~ 4 is 
better fit with a ~ —0.8, a value that was also found in the similar sized sample of 2; > 3 
quasars (Schneider, Schmidt & Gunn 1991). A change of continuum slope by 0.4 will result 
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in a 0.2 mag difference in the K-correction, or a 30% difference in quasar density. However, 
when we compare the observed results at 2; ~ 4 with z ~ 2 the situation is actually better 
than this calculated difference would imply. The low-redshift quasar surveys (such as 2dF) 
are selected based on the observed i?-band magnitude {Xeff ~ 4400A). At z ~ 2, the 
5-band corresponds to 1470 A in the rest-frame, almost the same as the case in which we 
select z ~ 4 quasars based on their i' band measurements {Xeff ~ 7600A). Therefore, both 
results are based on the rest-frame UV flux; the comparison is not strongly affected by the 
K-correction, as long as the two results are transferred to the rest-frame B-band in the 
same way. 

The presence of dust might bias the results on quasar evolution. Dust could arise along 
the line of sight in the damped Lya systems, or in the quasar environment, [ball fc Pel 



1(1993)1 studied the effect of dust in damped Lya systems on the determination of the quasar 
luminosity function. They found that 10% - 70% of quasars at z ~ 3 could be missed from 
optical surveys whose selection was based on B band magnitudes. This fraction cannot 
be calculated precisely, due to uncertainties in the dust content and dust distribution in 
damped Lya systems. But the effect of intervening dust must increase quickly with redshift, 
both because the total number of damped systems increases, and because the extinction 
curve rises sharply towards the rest-frame UV. The number density of quasars at z > 4 
could be significantly affected by dust along the line of sight, but it is not clear whether it 
is enough to change the general trend of the decline. 



High-redshift quasars are known to reside in dusty environments ( Pmont et al. 1996 



[Uarilli et al. 20001) , and the new millimeter and submillimeter surveys have revealed a 



population of high-redshift dusty galaxies (e.g. Parger et al. 1998]) . The existence of dust 



extinction changes the bolometric correction from the rest-frame UV luminosity observed 
in the optical surveys. If the bolometric correction is a strong function of redshift, the 
evolution of the bolometric luminosity function of quasars could be very different from that 
of optically selected quasars. However, we've seen that the a distribution of 2; ~ 4 quasars 
is quite similar to that of lower-redshift quasars, implying that the objects in our sample 
are no more reddened than those at low redshift. Thus in order to invoke dust to explain 
the decline of optical quasars with redshift, one would require a completely separate, very 
extincted population of objects whose numbers increase with redshift. 



Radio and X-ray selected samples are not strongly affected by dust extinction. [Shaver 



et al. (1996)1 found that the space density of high-redshift radio quasars also declines rapidly 



at 2 > 3. On the other hand, [IVliyaji, Hasinger fc Schmidt (200(T)[ found that there is no 



evidence for a decrease in the number density of RO SAT-selected X-ray quasars at 2; > 3, 
although the statistical significance of this result is only marginal. One difficulty is how to 
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normalize the results from different wavelength samples. At z ~ 3, the spatial density of the 
quasars in the [Miyaji, Hasinger fc Schmidt (200n]] sample is considerably higher than in the 



SSG sample, which implies that these X-ray quasars are intrinsically fainter objects than 
those in SSG. The redshift evolution is expected to depend on luminosity, and there are 
theoretical reasons to believe that the evolution is flatter at low luminosities (see below). 

Studies of the quasar luminosity function put important constraints on models of 
quasar evolution. Haiman & Loeb (1998) show that the observed quasar luminosity function 
can be fit with a theoretical model based on the Press-Schechter (1974) approximation, with 
reasonable assumptions about the relation between quasar luminosity and dark matter halo 
mass, and the duty cycle of quasar activity. Haiman & Hui (2000) and Martini & Weinberg 
(2000) further demonstrate that the duty cycle can be better constrained by fitting both 
the luminosity function and two-point correlation function of high-redshift quasars. Faint 
quasars reside in less massive dark matter halos in this model. They represent less rare 
peaks in the density field, and are expected to evolve more slowly than are the bright 
quasars. The luminosity function in this paper covers only two magnitudes in luminosity, 
which is not enough to make a detailed comparison with these models. A deep high-redshift 
quasar survey is needed in order to constrain the faint end slope of the quasar luminosity 
function and to compare with theoretical models. 

Observations of a 2; = 5.8 quasar ( Pan et al. 2000b|) demonstrates that the universe 
is already highly ionized at 2; ~ 6, by either the UV ionizing radiation from quasars or 
from star-forming galaxies. Assuming that the high-redshift quasar luminosity function has 
the same shape as at low redshift, Madau et al. (1999) estimated that the UV ionizing 
photons from high-redshift quasars is not sufficient to keep the universe ionized at z ~ 4. 
In this paper, we have shown that the shape of the high-redshift quasar luminosity function 
is shallower than at low redshift for bright quasars, implying even fewer ionizing photons 
from the these quasars than assumed in Madau et al. (1999). However, the faint quasars 
might contribute more UV ionizing photons than do bright quasars at high redshift. We 
currently have no knowledge of the faint end slope of the quasar luminosity function at 
z > 3. Therefore, a deep quasar survey is also needed in order to determine the nature of 
high-redshift ionizing background. 

The SDSS southern survey (see Paper HI) will image the Fall Equatorial Stripe (the 
same area of the sky as the sample used in this paper) 35 - 40 times in the five-year survey 
period. By co-adding these data, we will be able to select high-redshift quasar candidates 
down to i' ~ 22.5 at 2; ~ 4, and even fainter for low-redshift quasars. In Paper HI, we 
present 18 quasars down to ~ 21 selected from areas that have been observed twice during 
the SDSS commissioning. The high-redshift quasar luminosity function from the southern 
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survey will reach Mb ~ —24, close to the traditional boundary between quasars and Seyfert 
galaxies (Mg = —23). This sample will give a more complete description of the statistical 
properties of high-redshift quasars. 

The Sloan Digital Sky Survey (SDSS) is a joint project of the University of Chicago, 
Fermilab, the Institute for Advanced Study, the Japan Participation Group, The Johns 
Hopkins University, the Max-Planck-Institute for Astronomy, Princeton University, the 
United States Naval Observatory, and the University of Washington. Apache Point 
Observatory, site of the SDSS, is operated by the Astrophysical Research Consortium. 
Funding for the project has been provided by the Alfred P. Sloan Foundation, the SDSS 
member institutions, the National Aeronautics and Space Administration, the National 
Science Foundation, the U.S. Department of Energy, and Monbusho, Japan. The SDSS Web 
site is |http : //www . sdss . org/| . The Hobby-Eberly Telescope (HET) is a joint project of 
the University of Texas at Austin, the Pennsylvania State University, Stanford University, 
Ludwig-Maximillians-Universitat Miinchen, and Georg-August-Universitat Gottingen. The 
HET is named in honor of its principal benefactors, William P. Hobby and Robert E. Eberly. 
XF and MAS acknowledge support from Research Corporation, NSF grant AST96-16901, 
the Princeton University Research Board, and a Porter O. Jacobus Fellowship. DPS 
acknowledges support from NSF grant AST99-00703. 

A. Including the Selection Function in Lynden-Bell's C Estimator 

Lynden-Bell's (1971) C~ estimator calculates the marginal distribution of quasars in 
luminosity or redshift by constructing a comparable set, the luminosity and volume limited 
subsample for each object in the sample. The comparable set in luminosity is defined as: 



where Lj is the limiting luminosity at redshift Zj. For a survey with a complicated 
selection function such as ours, the limiting luminosity is not well-defined: the selection 
probability goes to zero smoothly at the survey limit; elsewhere, the selection probability 
is a complicated function of z, L and the spectral energy distribution of the object. In this 
case, the comparable set has to be redefined to include the contribution from this selection 
function. We follow the derivation of Petrosian (1992) for the C~ estimator, now including 
selection probabilities. 

Here we assume that the distribution of L and z are not correlated, thus the luminosity 
function can be written as: 



Ji — {j : Lj > Li, Lj < L j } , 



(Al) 



^iL,z)=ijiL)p{z). 



(A2) 
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The cumulative luminosity function is defined as $(L, z) — "^{L', z)dL', 

p{L,z) 



n{L, z)dL 



ir^^p{L,z)dL' 

where n{L, z) is the observed distribution and p(L, z) is the selection function. 



(A3) 



In order to get the marginal cumulative luminosity function, = ■ilj{L')dL' , 

we need to integrate over z. Since the quantity d(f){L) / (f){L) = d^{L, z)/^{L, z) does not 
depend on it is easy to show that 

#(^) ^ /o°" dz n{L, z)dL 

m J,-dzjrdL'^p{L,zy 

Integrating it gives: 

^/rN A i f°° S^dzn{L,z)dL \ 

ML) = ^exp ^ / dL '-^ ^^.7^^ ) , (A5) 

lo-^dzj-dL'^pmz)!' 

where ^4 is a constant. 

In the discrete case, the observed distribution is the sum of a series of 6 functions over 
individual objects: 

i 

Let us assume that the objects in the sample are sorted by their luminosities: 
Li > L2 > ... > Li-i > Li > ... > Ltv- Note that for the discrete case, is the step 

function: (piLj — e) = (j){Lj_i), and (piLj + e) = (f){Lj), when e <^ 1. Thus, one re- writes the 
integral from Lj — e to Lj + e, and regards p{L, z) as a constant over the integrated range: 



6{L - Lj)dL 



cl>{L,) = 0(L,_Oexpy , (A7) 

where Q{L) is a step function. We can define the comparable set for object j as: 

Jj ^{i:Li> Lj}, (A8) 

and 
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Note that de{x) = 5{x). We have: 

= ^(VO^Mll (A„) 

Thus, we have the C~ estimator for the cumulative luminosity function: 

0(L,) = 0(Li)n(l + l/Ar,) (A12) 

1=2 

for the case of an arbitrary selection function. 
Similarly, if we define: 

p(L- z) 

= V ^;7' ^ , where the sum extends over Zj < Zi, (A13) 
■ p{Lj,Zj) 

and the sample is sorted by redshift, Zi < ... < Zj^i < Zj < ... < Zn, the cumulative redshift 
distribution can be estimated as: 

k 

a{z,) = a{zi)f[{l + l/M,). (A14) 
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Table 1. Maximum Likelihood Solutions 



Parameter 


n = 1 


fl = 0.35 A = 0.65 




h = Q.5 


h = 0.65 


Differential, Eq. (26) 


*o (Mpc~^ mag~^) 




7.2lt^ X 10-« 


a 


0.76 ±0.29 


0.75 ± 0.29 




-2.58 ±0.23 


-2.58 ±0.23 


Cumulative, Eq. (25) 


A 


-7.24 ±0.19 


-7.31 ±0.19 


B 


0.48 ±0.15 


0.47 ±0.15 


C 


0.63 ±0.10 


0.63 ±0.10 
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Figure 1. The high-redshift quasar luminosity function derived from the l/K estimator, (a) 
is the result for the 0,-1 model, corrected for the average selection function; (b) is the 
result for the Q — 1 model, corrected for the selection function using each quasar's SED 
type; (c) and (d) are for the A-model with the two methods of selection function correction. 
The dashed lines are the maximum likelihood solutions. 
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Figure 2. The high-redshift quasar luminosity function derived using Lynden-Bell's (1971) 
C~ estimator, (a) and (c) show the marginal differential luminosity distribution as a function 
of magnitude, for both the — 1 and the A-model; (b) and (d) show the marginal redshift 
evolution: the spatial density of quasars at M1450 < —25.5 as a function of redshift. The 
dashed lines are the maximum likelihood solutions. 
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Figure 3. The evolution of the quasar spatial density at Mb < —26 is compared with previous 
studies. The SDSS points shown are the results using the C~ estimator, and the dashed line 
is the maximum likelihood solution. The SSG points shown are the spatial density calculated 
with the 1/Va method. The low-redshift result is the best-fit model from the 2dF survey 
(Boyle et al. 2000). 
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Figure 4. The cumulative luminosity function at z ~ 4.3, compared with previous studies. 
The SDSS result shown (shaded area) is based on the 1/K estimator for A.O < z < 4.5. The 
SSG result is their best-fit model (Eq. 35). For the z — 2 luminosity function, we use the 
best- fit model from the 2dF survey (Boyle et al. 2000). 
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Figure 5. The slope of the high-redshift quasar luminosity function at the bright end. The 
SDSS result gives i/^iL) oc L~^-^^ (solid line), considerably shallower than the low-redshift 
luminosity function ■0(L) oc L~^-^^. The dashed line in the plot is the maximum likelihood 
solution from the SDSS sample if the slope of the luminosity function is forced to have 
P = -3.43. 



